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First-principles study on field evaporation for silicon atom on Si(OOl) surface 
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The simulations of field-evaporation processes for silicon atoms on various Si(001) surfaces are 
implemented using the first-principles calculations based on the real-space finite-difference method. 
We find that the atoms which locate on atomically flat Si(00f ) surfaces and at step edges are easily 
removed by applying external electric field, and the threshold value of the external electric field for 
evaporation of atoms on atomically flat Si(00f) surfaces, which is predicted between 3.0 and 3.5 
V/A, is in agreement with the experimental data of 3.8 V/A. In this situation, the local field around 
an evaporating atom does not play a crucial role. This result is instead interpreted in terms of the 
bond strength between an evaporating atom and surface. 

PACS numbers: 



I. INTRODUCTION 

In recent years, scanning tunneling microscopy exper- 
iments in the atomic-scale manipulation of material sur- 
faces have attracted great attention due to its potential 
to create artificial surface nanostructures. It is now pos- 
sible to remove adsorbed atoms and molecules from sur- 
faces or to deposit them on surfaces by applying voltage 
pulses Field evaporation and deposition are ther- 
mally activated processes in which rate constants can be 
parameterized according to the Arrhenius formula [3j . In 
addition, the electric field between probing tip and sur- 
face plays a crucial role during these processes, since the 
activation energy for field evaporation is considered to de- 
pend on this electric field. In the past, phenomenological 
models have been used to evaluate the activation energy 
for field evaporation of surface atoms Q- However, in 
order to understand this phenomenon properly from a 
microscopic point of view, it is mandatory to calculate 
self-consistently both the electronic charge distribution 
and the resulting electrostatic field. 

In this paper, we carry out first-principles molecular- 
dynamics simulations to explore the activation energies 
and threshold values of external electric fields for the 
evaporation of surface atoms from Si(001) surfaces. Five 
years ago, Kawai et al. 3 determined the adiabatic po- 
tential curves for adsorbed atoms on the atomically flat 
Si(001) surface, however, the maximum external electric 
field of ~ 2.6 V/A applied in their simulation is rather 
small compared to that in experiments using field ion 
microscopy (3.8 V/A) 0. We now determine activation 
energies for silicon atoms on Si(001) surfaces under the 
stronger electric fields of 3.0, 3.5, and 4.0 V/A. More- 
over, we examine differences in activation energies and 
threshold values of external electric fields among silicon 
atoms evaporating from various Si(001) surfaces, e.g., the 
atomically flat surface, the surface having a step, and the 



atomically flat surface with an adsorbed silicon atom. We 
find that the threshold value for field evaporation of the 
atom on the atomically flat Si(001) surface, which is be- 
tween 3.0 and 3.5 V/A, is the lowest and that for the 
surface atom in the atomically flat surface, which is over 
4.0 V/A, is the highest. These results will be explained 
in terms of binding energies between evaporating atoms 
and surfaces. 

The rest of the paper is organized as follows: in Scc.lTTl 
we describe briefly the computational details of our cal- 
culations. The results are presented and discussed in 
Sec. IIIII We conclude our findings in Sec. IIVI 



II. COMPUTATIONAL DETAILS 
A. Method 

Our first-principles molecular-dynamics simulations 
are based on the real-space finite-difference method [j| 
with incorporation of the timesaving double-grid tech- 
nique 6|. Compared with the plane- wave approach, the 
real-space finite-difference method is much simpler to im- 
plement while keeping a high degree of accuracy. More- 
over, the real-space calculations eliminate the serious 
drawbacks of the conventional plane- wave approach such 
as its inability to describe strictly nonperiodic systems: 
in the case of simulations under external electric fields, 
the periodic boundary condition gives rise to saw-tooth 
potential, which sometimes leads to numerical instability 
during the self-consistent iteration, while the real-space 
finite-difference method can exactly determine the po- 
tential by the external electric field as a boundary con- 
dition and is free from involving the saw-tooth potential. 
We adopt the nine-point finite-difference formula for the 
derivative arising from the kinetic-energy operator of the 
Kohn-Sham equation in the density functional theory Q. 
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TABLE I: Total energies as the atom is displaced at sites A-G. 
The zero of total energy is chosen to be that at site E. 



bite 


Total energy (eV) 


A 


1.42 


B 


1.71 


C 


0.98 


D 


1.51 


E 


0.00 


F 


1.14 


G 


1.28 



We take a cutoff energy of 23 Ry, which corresponds to 
a grid spacing of 0.65 a.u., and a higher cutoff energy of 
210 Ry in the vicinity of nuclei with the augmentation of 
double-grid points |6J. The norm conserving pseudopo- 
tential [8|, la is employed in a Kleinman-Bylander non- 
local form [lOj . Exchange-correlation effects are treated 
with local density approximation [llj . 



B. Models 



Figure ^ shows the top views of the Si(001) surfaces 
calculated here. We employ a technique that involves 
the use of a supercell whose size is chosen as L x =43.52 
a.u., Lj,=14.51 a.u. and L z =32.64 a.u., where L x , L y 
and L z are the lengths of the supercell in the x, y and z 
directions, respectively. Here the direction perpendicu- 
lar to the surface was chosen as the z direction. In order 
to eliminate completely unfavorable effects of atoms in 
neighbor cells which are artificially repeated in the case 
of the periodic boundary condition, we impose the non- 
periodic boundary condition of vanishing wave functions 
out of the supercell in the z direction, while we adopt the 
periodic boundary condition in the x and y directions. 
The supercell contains five silicon layers, the lowest of 
which is terminated by hydrogen atoms (i.e., the thin 
film model). The atoms in the three topmost silicon lay- 
ers are fully optimized by the first-principles molecular- 
dynamics simulation. In the cases of an atomically flat 
surface [terrace model in Fig. ^ (a)] and a surface having 
a step [Sa step model in Fig. HJ(b)], we employ previ- 
ously reported atomic configurations 0. On the other 
hand, for the atomically flat surface where an atom is 
adsorbed [adatom models in Figs. Q]( c ) an d (d)], we de- 
termine the minimum-energy atomic configuration. The 
atom is initially placed at each of the sites of A, B, C, 
D, E, F, and G in Fig. ^with an appropriate distance 
from the surface, and then the forces on the atoms are 
relieved during the structural optimization. The total 
energies are collected in Table H] The most stable site is 
found to be site E in Fig. ^ (d) , and hereafter we use this 
configuration as the adatom model. 



FIG. 1: Top views of the three topmost layers of Si(001). (a) 
the atomically flat surface (terrace model), (b) the surface 
having a Sa step (Sa step model), (c) and (d) the surfaces 
where an atom is adsorbed (adatom models). Closed circles, 
opened circles with solid curves, and opened circles with dot- 
ted curves represent atoms on the top, second, and third lay- 
ers, respectively. Large (small) circles represent the upper 
(lower) atoms in buckled dimers. 



III. RESULTS AND DISCUSSION 

A. Field evaporation of surface atoms 

We first apply an external electric field F of 3.0 V/A 
along the z direction and find that no surface atoms evap- 
orate in any models. When the external electric field 
F is increased to 3.5 V/A, the adsorbed atom in the 
adatom model evaporates [see Fig. El (a)]. On the basis 
of these results, it is evident that the threshold value of 
the external electric field for evaporation is 3.0-3.5 V/A. 
This value is in agreement with the experimental data of 
3.8 V/A obtained using a field ion microscope. In the 
case of ^=4.0 V/A, the two upper atoms of the buckled 
dimers located at the step edges simultaneously evapo- 
rate in the S A step model [see Fig. 0(b)] 13]. The 
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surface atoms in the Sa step model are preferentially 
removed from the step edges. No surface atoms in the 
terrace model evaporate at all below F—A.O V/A. 



TABLE II: Activation energies (ineV) Q. The values at zero 
field are the binding energies between the lifted-up atoms and 
the surfaces. 



terrace 
Sa step 
adatom 



0.0(V/A) 3.0(V/A) 3.5(V/A) 4.0(V/A~ 



6.03 
5.86 
5.04 



1.45 
1.05 
0.09 



0.52 
0.13 
0.00 



0.24 
0.00 
0.00 



FIG. 2: Isosurfaces of the charge density for (a) the adatom 
model under F=3.5 V/A and (b) the Sa step model under 
F=A.O V/A. Light and dark balls represent the surface atom 
and the atoms in the silicon surface, respectively. 



B. Activation energy 



One of the potentially important applications of the 
field-evaporation process is the direct determination of 
the binding strengths of the surface atoms from the ex- 
ternal electric field required for their removal. We now 
evaluate the activation energies of surface atoms for field 
evaporation during the lifting up of the surface atoms. 
Table ITT1 shows the activation energies for various surface 
models The activation energies in the absence of 

the external electric field, which are equal to the binding 
energies between evaporating atoms and surfaces, are ob- 
served to be 6.03, 5.86, and 5.04 eV for the terrace, Sa 
step, and adatom model, respectively. These energies be- 
come lower as the external electric fields are increased. 
The threshold value of the external electric field varies 
with the bond strength of the surface atom, and the rate 
constants of field evaporation depend on the activation 
energies according to the Arrhenius formula [3j; surface 
atoms in both the Sa step model and the adatom model 
are easily removed compared with the surface atom in 
the terrace model. This situation may lead to the devel- 
opment of technology to create atomically flat surfaces 
by applying the external electric fields to surfaces and 
adjusting the temperature. 



FIG. 3: Electronic density shift p(r, F)—p(r, F = 0) at F=3.5 
V/A represented on the cross section of the ribbon con- 
taining the thick line A-B. The contour spacing is 6.7 elec- 
tron/supercell. Solid (dotted) curves represent nonnegative 
(negative) values. The large and small balls indicate the 
atomic positions on and above the cross section, respectively. 



C. Local-field enhancement 

Figure [21 shows the electron density shift p(r, F) — 
p(r, F — 0) due to the application of the external elec- 
tric field of -F=3.5 V/A[l5|]. The overall charge around 
the surface atoms decreases, and the atom is expected 
to be a positive ion when it evaporates. We then cal- 
culate the external electrostatic field as the difference 
between the total electrostatic field in the presence of 
an external electric field and that in the absence of it, 
d(V e ff(r, F) — V e ff(r, F — 0))/dz, where V e // is the sum 
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difference in external electrostatic field around the evap- 
orating atom is observed. Therefore, we cannot conclude 
that the local field around the evaporating atom plays 
a crucial role in the situation mentioned above in Sec. 

MM 



IV. CONCLUSION 



FIG. 4: Difference in external electrostatic field between 
F=0.Q V/A and F=3.5 V/A. The contour spacing is 0.5 
V/A. The meanings of the symbols are the same as those in 
Fig. H 

of the external, Hartree, and exchange-correlation po- 
tential functions. Figure 21 shows the counter plots of 
the difference in external electrostatic field 0] . The ex- 
pulsion of the external electronic field from inside the 
surface can be clearly recognized in all models. How- 
ever the local-field enhancement occurs not around the 
evaporating atom but above it, and no further significant 



We have studied the field-evaporation process of sur- 
face atoms by using various surface models. The thresh- 
old value of the external electric field for evaporation of 
the surface atom on the atomically flat Si(001) surface, 
which is between 3.0 and 3.5 V/A and is in agreement 
with the experimental result of 3.8 V/A, is the lowest 
among the models in this study, and no surface atoms 
of the atomically flat Si(001) surface evaporate below 
the external electric field of 4.0 V/A. When the exter- 
nal electric field is applied to the surface having a Sa 
step, the atoms which are located at the step edges are 
preferentially removed. The threshold value of the exter- 
nal electric field obviously depends on the binding energy 
between evaporating atoms and surfaces. Moreover, as 
the external electric field increases, the activation energy 
for field evaporation becomes lower. On the basis of the 
discussions above, we conclude that the field-evaporation 
process can be applied to surface-flattening techni ques 
|l6j or atomic-scale manipulation methods 0, ITsl Il9f 
by adjusting the strength of the external electric field 
and the temperature. 



Acknowledgements 

This research was partially supported by the Ministry 
of Education, Culture, Sports, Science and Technology, 
Grant-in-Aid for COE Research (Grant No. 08CE2004) 
and Young Scientists (B) (Grant No. 14750022). The 
numerical calculation was carried out by the computer 
facilities at the Institute for Solid State Physics at the 
University of Tokyo, and Okazaki National Institute. 



[1] See for example D.M. Eigler and E.K. Schweizer in Na- 
ture 344 524 (1990), and I.-W. Lyo and P. Avouris in 
Science 253 173 (1991). 

[2] T.T. Tsong, Atom Probe Field Ion Microscopy, (Cam- 
bridge University Press, Cambridge, 1990), and refer- 
ences are therein. 

[3] S. Arrhenius, Z. Phys. Chem. 4, 226 (1889). 

[4] K. Watanabe and T. Satoh, Surf. Sci. 287-288, 502 
(1993); T. Kawai, and K. Watanabe, Surf. Sci. 357-358, 
830 (1996); T. Kawai, and K. Watanabe, Surf. Sci. 382, 
320 (1997); T. Kawai, K. Watanabe, and K. Kobayashi, 
Ultramicroscopy 73, 205 (1998). 

[5] J. R. Chelikowsky, N. Troullier, and Y. Saad, Phys. Rev. 



Lett. 72, 1240 (1994); J. R. Chelikowsky, N. Troullier, 
K. Wu, and Y. Saad, Phys. Rev. B 50, 11355 (1994). 

[6] T. Ono and K. Hirose, Phys. Rev. Lett. 82, 5016 (1999). 

[7] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 

(1964) ; W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 

(1965) . 

[8] N. Troullier and J.L. Martins, Phys. Rev. B 43, 1993 
(1991). 

[9] We used the norm-conserving pseudopotentials NCPS97 
constructed by K. Kobayashi. See K. Kobayashi, Corn- 
put. Mater. Sci. 14, 72 (1999). 
[10] L. Kleinman and D.M. Bylander, Phys. Rev. Lett. 48, 
1425 (1982). 



5 



[11] J. P. Perdew and A. Zunger: Phys. Rev. B 23, 5048 
(1981). 

[12] We employed the atomic geometries of the Si(001) sur- 
face reported by A. Ramstad, G. Brocks, and P.J. Kelly 
in Phys. Rev. B 51, 14504 (1995) and D.J. Chadi, Phys. 
Rev. Lett. 59, 1691 (1987) as initial atomic configura- 
tions for the terrace and Sa step models, respectively, 
and then carried out structural optimizations. 

[13] There are two Sa step edges in the supercell, because 
we imposed periodic boundary condition in the direc- 
tions parallel to the surface. Each of two upper atoms of 
the dimers is located slightly different situation, however 
deferences in activation energies for the field evaporation 
between them are negligibly small. 



[14] For the Sa step model, we depict the activation energies 
of the center atom in Fig.|H(b). 

[15] In Figs. [3] and |1J the charge distribution p(r, F = 
3.5V/A) and potential V eff {r,F = 3.5V/A) are gener- 
ated using the stable atomic configuration at F=Q.O V/A. 

[16] K. Sudoh and H. Iwasaki, Jpn. J. Appl. Phys. 39, 4621 
(2000). 

[17] D.M. Eigler and E.K. Scheweizer, Nature 344, 524 
(1990). 

[18] I.-W. Lyo and P. Avouris, Science 253, 173 (1991). 
[19] D.M. Eigler, CP. Lutz and W.E. Rudge, Nature 352, 
600 (1991). 



This figure "figl.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0307683vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0307683vl 



This figure "fig3.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0307683vl 



This figure "fig4.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0307683vl 



